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Motivated by the recent experimental work (A. Widera, et al, Phys. Rev. Lett. 95, 19045), 
we study the collisional spin dynamics of two spin-2 87 Rb atoms confined in a deep optical lattice. 
When the system is initialized as 0, 0), three different two-particle Zeeman states are involved in 
the time evolution due to the conservation of magnetization. For a large magnetic field B > 0.8 
Guass, the spin coherent dynamics reduces to a Rabi-like oscillation between the states jO, 0) and 

1, —1). However, under a small magnetic field, a general three-level coherent oscillation displays. 
In particular, around a critical magnetic field B c ~ 0.48 Guass, the probability in the Zeeman states 

2, —2) exhibits a novel quantum beat phenomenon, ready to be confirmed in future experiments. 
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Spinor Bose-Einstein condensates (BEC) in purely op- 
tical traps has opened a new direction in study of con- 
fined dilute atomic gases 0,0,|il,GllIlEI0! an( A many fas- 
cinating phenomena originatingfrom the spin degrees of 
freedom have been observed |sL Ifl ITfl ITU . Among 

them, coherent spin-exchange dynamics induced by inter- 
atomic collisions has been investigated in several recent 
experiments with BEC condensates in an optical lattice 
IHJ LLjJ LUi LLJ] j where coherent control of the evolution 
with a magnetic field to apply different phase shifts to 
the spin states has also been demonstrated fl5| . 

More recently, coherent spin-mixing oscillations in a 
Mott insulating state of hyperfine spin-2 87 Rb atoms in 
deep optical lattices have been reported in a system 
of many isolated of a pair of atoms on each lattice site. A 
weakly damped Rabi-like oscillation between two-particle 
Zeeman states with equal magnetization displays in the 
large magnetic field, the oscillation frequency and am- 
plitudes are precisely determined by the spin-exchange 
couplings and the second order Zeeman shifts. However, 
when the external magnetic field is weak, three different 
two-particle Zeeman states with equal magnetization are 
involved, and a general feature of three-level coherent os- 
cillations are expected among the possible spin collisional 
processes. 

In this paper, we present the detailed calculations of 
the coherent evolutions of a pair of spin-2 87 Rb atoms 
in a deep optical trap. When the system is initialized 
as |0,0), for a large magnetic field typically B > 0.8G, 
the spin coherent dynamics reduces to a Rabi-like oscilla- 
tion between the states 1 0, 0} and |1, —1), which has been 
observed in the excellent experiment 16]. However, for 
a weak magnetic field, a three-level coherent oscillation 
displays. Around a critical magnetic field B c 0.48G 
for the typical parameters, it is more interesting that the 
probability in the Zeeman state |2,— 2) exhibits a novel 
quantum beat phenomenon, which is possible to be ob- 
served in further experiments. 

For a pair of spin-2 boson atoms in a deep optical trap, 



the general interaction is given by 

V = ffoPo +52P2 +34P4, (1) 

where gp = 4Trh 2 ap/M J dr|</?| 4 , with M the mass of 
the atom, a F the s-wave scattering lengths in the total 
spin-F channel, Pf the projection operator for a total 
hyperfine spin-F state with Po + P2 + P4 = 1, and ip is 
the spin-independent spatial wave function of the ground 
state in the potential trap. 

For a pair of atoms, one could use the relation Si-S2 = 
— 6P0 — 3P2 + 4P4 to rewrite the interaction into 

K = co + c 1 SrS 2 + c 2 Po, (2) 

with c = (4.92 + 3p 4 )/7, c x = (g 4 - ch)/7, c 2 = 
(7g - 10g 2 + 2,g A )/l and P = 0+6/10 H U, where 
0t = 2b\b*_ 2 — < 2b\b^_ 1 + frJ&J creates a pair of atoms in 
the spin singlet when applied to the vacuum, and 6^ the 
creation operator. Since the total spin angular momen- 
tum and its z-component are conserved, the interaction 
can be diagonalized with the eigenvalues Ej = ± = 10ci, 

Ej-2 = 3ci, and -Ej=o = C2, where j is the good quan- 
tum number of the total angular momentum. The corre- 
sponding eigenstates are described by |^> j m ). It should 
bc noted that one could not directly obtain the atom 
populations in these eigenstates experimentally. Instead, 
the populations of atom pairs in Zeeman states 1 7711,7712) 
could be easily obtained by absorption imaging after a 
few milliseconds of time-of-flight (TOF) with a magnetic 
gradient field. 

When the atom pair is initialized as the Zeeman state 
1 0, 0) , other two-particle Zeeman states in the subspace 
of mi + = are involved in the time evolutions due 
to the conservation of the total spin z-component mag- 
netization. Introducing the following notation 
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we can explicitly express the two-particle Zeeman states where the coefficients in front of the eigenstates 



in terms of linear combinations of the eigenstates |0 ■ m ). 
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^|04,O> + Y^|^2,O> + y||0O,O>. ( 4 ) 
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for each Zeeman states |mi,m2) are just given by the 
Clebsh-Gordon (C-G) coefficients. Therefore, the wave 
function of the system at any time t can be expressed 



\i/> t ) = C„(t) |0, 0) + d(t) |1, -1) + C 2 (t) |2, -2) , 



with 
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The corresponding probabilities in each Zeeman state are 
easily derived 

473 72 4 36 

Pn — 1 COSCthtH COSCJt^H COSUJ-it, 

1225 245 35 175 3 ' 

436 48 4 48 
-Pi = cosw-iiH cosujnt coswai, 

1 1225 245 35 175 3 ' 

316 24 8 12 

Pn — COS Wit COS LUnt H COSLU-lt, 

2 1225 245 35 175 3 ' 

where w\ — E 4 , — E 2 fi, uj 2 = E 2 , — E , , and w 3 = 
E^fi — £uo- Clearly, the spin dynamics of the system 
displays a three-level coherent oscillation with three dif- 
ferent frequencies and comparable amplitudes. We have 
plotted these probabilities in Fig.l. We would like to em- 
phasize that the oscillation frequencies are determined 
by the eigenvalues, instead of the off-diagonal matrix el- 
ements, of the 3x3 model Hamiltonian matrix. 

In the presence of a magnetic field, the linear and 
quadratic Zeeman terms should be included into the in- 
teraction Hamiltonian [13 



H = V - p(S lz + S; 



1z ) 



q(S 



(6) 



where p = fi B B/2h, the ground state hyperfine split- 
ting uififs ~ 27r x Q.835GHz for 87 Rb atoms, and q = 
p 2 /uj h f s rj 450B 2 Hz/G 2 . In this case, the total spin 
angular momentum is no longer conserved, but its z- 
component is still conserved. Therefore, we can work in 
the subspace of mi + 7712 = 0, in which the second or- 
der Zeeman shift can influence the spin dynamics. In the 
complete bases of |0, 0), |1, —1), and |2, —2), the effective 
Hamiltonian could be expressed as: 
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FIG. 1: Probabilities of the two-particle Zeeman states in the 
absence of magnetic field. 



where we have used ao = 89.4as, a 2 = 94.5as, and 04 
106.0a B for 87 Rb atoms 0, and J - 1 ' i! 



M 



■jdrM' 

is typically chosen as 129.02ffz. To solve this Hamilto- 
nian, one has to employ the numerical method. Using 
the harmonic wave function to approximate ip, one could 
estimate 



Ja B 



ash 2 



irM 
4:.08v 3/4 Hz 



-(^V^M) 3/2 
7t n 



(8) 



where Vq = is the depth of the lattice potential. 

In the Mott insulating regime, v is large and we simply 
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FIG. 2: Eigenvalues of the effective model in (a) and the 
corresponding three oscillation frequencies in (b) as a function 
of magnetic field. 



choose v = 100. In Fig. 2a, the eigenvalues of the effective 
Hamiltonian (E\ > E 2 > E3) have been presented as a 
function of the magnetic field. 

In the presence of the second order Zeeman shift, three 



two-particle Zeeman states can still be expressed in terms 
of three eigenstates 
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where nine coefficients are calculated numerically and are 
plotted as a function of the magnetic field in Fig. 3a. It 
can be seen that Uxi and C/33 become very small and re- 
duce quickly to zero as the magnetic field is increased. 
When U\\ = 0, the initial state is just a superposition 

of two eigenstates only |0,0) ~ U12 



U 
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it will not jump into third eigenstate Eij. Similar to 

the previous procedure, the probabilities of the Zeeman 
states are obtained as the form of 

P = A 00 + A i cos wit + A 2 cos uj 2 t + A 03 cos cu 3 t, 
Pi = A10 + An cos wit + A12 cosu> 2 t + A13 coso^i, 
P2 = A 2 q + A21 cosujit + A22 coso^ + A23 cosLu 3 t, 

with the amplitudes as combinations of the coefficients 
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and the frequencies are given by 

loi = Ei — E2, L02 = E2 — E3, LU3 = Ei — E3. 

There are simple relations between the amplitudes of 
each oscillations and the frequencies: 

Aon + Ai n + A 2n = S n>0 , w 3 =wi+w 2 - (11) 

As the magnetic field B gradually increases, the fre- 
quencies oji of the spin oscillations are determined by 
the eigenvalues E n , while the amplitudes of each fre- 
quency oscillation are determined by the coefficients U%j . 
By focusing on the amplitudes A mn for each oscillation 
frequency as the functions of the magnetic field shown 
in Fig. 3b, we can clearly demonstrate how the system 
crossovers from a three-level (|2,— 2), |1,— 1), |0, 0)) co- 
herent oscillation without any selection rules into a two- 
level (|1,— 1), 1 , 0) ) Rabi-like oscillation. 



When the large magnetic field is typically larger than 
0.8G, which is a reasonable magnitude compared to the 
experimental parameters |l6j . all amplitudes A mn begin 
to fall quickly down to zero except for Aq2 and Ai 2 . 
Thus, P 2 (t) could be neglected and Po and Pi are re- 
duced into an oscillation with a cosine form, i.e., the 
system oscillates between |0, 0) and |1,— 1) like a Rabi 
oscillation. In Fig. 4, we have plotted the probabilities of 
the three Zeeman states at B — 1.2G. To find out the 
more detailed reason, one must look into the coefficients 
Uij shown in Fig. 3b: Un and U33 become very small and 
reduce quickly to zero, A m i and A m 3 (m = 1,2,3) in- 
cluding Un fall down to zero, and ^22 including U 33 falls 
down to zero too. In the extreme limit of q/J S> 1, all 
Uij are zero except U\ 3 , U 22 , U31 equal to 1. Then the 
initial state |0, 0) just corresponds to one eigenstate, so 
there will be no spin-exchange collisional dynamics. 

However, in the weak magnetic field regime, the ampli- 
tudes of three frequencies in Po, Pi and P 2 are all finite 
and comparable, displaying the three-level coherent oscil- 
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experimental predictions are given in the weak magnetic 
field limit, and we have also demonstrated the presence 
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FIG. 3: Expansion coefficients of the two-particle Zeeman 
states in terms of eigenstates Uij in (a) and the amplitudes of 
each frequency oscillation Aij in (b) as a function of magnetic 
field. 



lations. As B is gradually increasing, the different eigen- 
values grow with different slopes. Around a critical mag- 
netic field B c fh 0.48G, ui and uj 2 are very close to each 
other, Aqx ~ A 02 ~ 0.15, A Q3 ~ 0.35, A n ~ A 12 ~ 0, 
A 13 ~ -0.42, A21 ~ ^22 -0.15, and A 23 ~ 0.07. 
Since A23 is small with compared to the beat ampli- 
tude (roughly 2A 2 i or 2^22), there appears a clear quan- 
tum beating phenomenon presence in the probability of 
the Zeeman state |2, — 2). Meanwhile, the correspond- 
ing probability Po(^) exhibits a combination of a beating 
oscillation and a sinusoidal one with the frequency 2u>, 
because the sinusoidal amplitude ^4o3 is comparable to 
the beat amplitude (roughly 2^4 i or 2A 02 ). Pi(t) just 
corresponds to a single frequency oscillation. In Fig. 5, 
we have plotted the three probabilities at B — 0.47G. 
Considering that this novel phenomenon is very sensitive 
to the external magnetic field, the critical point might be 
accurately determined by fine-tuning the external field 
experimentally. 

In conclusion, we have studied the collisional spin dy- 
namics of an isolated spin-2 87 Rb atom pairs confined in 
a deep optical lattice. Although our calculations could 
not consider the damping of the oscillations, some clear 
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FIG. 4: Probabilities of the two-particle Zeeman states in the 
large magnetic field B = 1.2G. 
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FIG. 5: Probabilities of the two-particle Zeeman states at 
B — 0.47G around the critical magnetic field. 



of a crossover from the three-level to two-level coherent 
spin dynamics. In particular, when the system is ini- 
tialized as 1 0, 0) and the magnetic field B c w 0.48G, the 
probability in the two-particle Zeeman state |2,— 2) will 
exhibit a quantum beat phenomenon, which is ready to 
be confirmed in future experiments. 
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